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We study the out-of-equilibrium large time dynamics of a gaussian polymer chain in a quenched 
random potential. The dynamics studied is a simple Langevin dynamics commonly referred to as the 
Rouse model. The equations for the two-time correlation and response function are derived within 
the gaussian variational approximation. In order to implement this approximation faithfully, we 
employ the supersymmetric representation of the Martin-Siggia-Rose dynamical action. For a short 
ranged correlated random potential the equations are solved analytically in the limit of large times 
using certain assumptions concerning the asymptotic behavior. Two possible dynamical behaviors 
are identified depending upon the time separation- a stationary regime and an aging regime. In the 
stationary regime time translation invariance holds and so is the fluctuation dissipation theorem. 
The aging regime which occurs for large time separations of the two-time correlation functions 
is characterized by history dependence and the breakdown of certain equilibrium relations. The 
large time limit of the equations yields equations among the order parameters that are similar to 
the equations obtained in the statics using replicas. In particular the aging solution corresponds 
to the broken replica solution. But there is a difference in one equation that leads to important 
consequences for the solution. The stationary regime corresponds to the motion of the polymer 
inside a local minimum of the random potential, whereas in the aging regime the polymer hops 
between different minima. As a byproduct we also solve exactly the dynamics of a chain in a 
random potential with quadratic correlations. 

PACS numbers: 05.40.-a, 05.70.Ln, 36.20.Ey, 75.10.Nr 

I. INTRODUCTION 

The behavior of polymer chains in random media attracted much interest in recent years because of its relevance 
and applications in diverse fields 0, ■ Besides elucidating the properties of the polymers chains themselves which is 
of much interest in physical chemistry |j| and biology 3 , this problem is directly related to the statistical mechanics of 
a quantumparticle in a random potential Q , the behavior of flux lines in superconductors in the presence of columnar 
defects 0,0) and the problem of diffusion in a random catalytic environment 0. It was found in Refs. [UElEl^SEi 
that a very long Gaussian chain, immersed in a random medium with very short range correlations of the disorder, 
will typically curl up in some small region of low potential energy. The polymer chain is said to be localized and for 
long chains the radius of gyration or the end-to-end distance becomes independent of chain length (R 2 ~ L°). 

Both heuristic arguments and a variational solution of the problem using replicas 01 > yielded the dependence 
of the size of the trapped polymer on the variance of the random potential (g), and the logarithm of the volume 
of the medium (V) such that R ~ (gin y) -1 /! 4 -^) f or 1 < d < 4. The breaking of replica symmetry was crucial 
to the derivation of the subtle InV dependence. In a related paper it was found that a quantum particle in a 
random environment exhibits glassy behavior at low temperatures. The low temperature limit for a quantum particle 
translates into the long chain limit for polymers. It implies that the free energy landscape of a long chain is typically 
very complicated and possesses many metastable states. In a recent publication [l3| we further utilized this mapping 
to give a physical interpretation of the localization and glassy behavior of a polymer in a random potential by making 
a connection with the localization of a quantum particle in a disordered medium of finite volume. Subsequently, 
we treated the case of random obstacles as opposed to a random potential, jl5| and finally included the effect of a 
self-avoiding interaction. 0, . 

Recently methods have been developed to study analytically the large time, non-equilibrium behavior of glassy 
systems. [3^3123 Directed polymers and manifolds have been investigated at the mean field level, and the solution 
exhibits two asymptotic time regimes: a stationary dynamic regimes at large but similar times and a slow aging 
regime for large and widely separated times. This large-time solution of the dynamical equations has many features in 
common with the replica-symmetry-breaking (RSB) solution of the corresponding equations of the statics, although 
replicas are not actually used and the limit n — > is avoided. But the equations, in particular for the case of 1-step 
RSB, are not all the same in the large time limit as the equilibrium equations. This leads to a situation that the 
system fails to reach the ultimate equilibrium state starting from an arbitrary initial condition, and ends up in a state 
with higher free energy than the one found at equilibrium. Thus the importance of the dynamical approach is twofold: 
to make contact with experiments that exhibit slow relaxation and aging effects |21j and also to serve in some cases 
as an alternative to the replica approach and the n — > limit, although this has to be taken with a grain of salt as 
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indicated above. 

It is our goal here to extend the previous treatment of large time dynamics for particles and directed manifolds in 
quenched random environment to the case of real polymer chains. The difference between the case of directed 

polymers and "real" polymers is mainly in the form of the random potential. If s denotes the bead number of the 
chain, for "real " chains different beads at the same spatial locations should be exposed to the same external potential, 
whereas for directed polymers different segments always feel a different random potential even at the same transverse 
position. This is easily made clear for flux lines in a superconductors. If we have point disorder different in different 
x-y planes when moving along the z-direction then the problem corresponds to directed polymers in a quenched 
random potential. On the other hand if we have randomly positioned columnar defects i.e. the random potential is 
independent of z, then when projected on a single plane we see that the system maps to "real" polymer chains in a 
quenched random potential. In addition a real polymer might have self-avoiding interactions among different beads, 
but these will not be considered in this paper. 

In this paper we consider the Langevin dynamics of a single gaussian chain embedded in a quenched random 
potential. This Langevin dynamics is referred to in the literature as the Rouse model for a polymer, and is the 
simplest dynamics 22] . For polymers in solutions a more realistic dynamics that reproduces more accurately the 
experimental results is the Zimm model that takes into account the effect hydrodynamic interactions effects. This 
kind of dynamics will not be considered in this paper and is a project for future research. Our goal here is to consider 
the simplest model that renders itself to an approximate analytical solution, of the l ong times d yna mics. Various 
treatments of polymer dynamics in a random potential have been considered before |23l 124. l25l |26j using various 
approximations, renormalization group treatment, and/or numerical solutions of approximate equations of motion. 
These papers implicitly assume time translation invariance (TTI) of the dynamical correlation function. They try 
to extract the behavior of the center of mass diffusion coefficient at short and large times, and their conclusions are 
not always in total agreement. Only recently the possibility of two-time dependence of the dynamical correlation 
function, i.e. aging phenomena has been explored numerically by Monte Carlo simulations, |27| but the dependence 
on the waiting time is not reported in detail in this paper. 

Our goal is also to make contact with the previous treatment of directed polymers in random potential and by 
using similar methods the difference between real and directed polymers will be elucidated. In addition, we derive 
the equations for a random potential with general correlations, either short or long range. The equations are derived 
using the so called gaussian variational approximation which is a kind of a mean field treatment. This approximation 
was first introduced by Fey nman in his studies of superfluid helium and in the context of random systems was first 
introduced in Refs. |3lLl34j . In order to implement it faithfully such as to preserve correctly all the symmetries of the 
dynamical equations we use the supersymmetric (SUSY) formulation of the dynamics p9j. implement the variational 
approximation, and then disentangle SUSY to produce coupled integro-differential e quat ions for the correlation and 
response functions. Here we follow steps similar to those used by Cugliandolo et al. |l8t Il9| and Konkoli et al. [2(j 
in their treatment of directed manifolds and random heteropolymers respectively. But because the random potential 
is implemented differently in our case, i.e. different beads at the same spatial position feel the same potential, the 
resulting equations are different. 

In order to solve the equations analytically at large times we make some assumptions about the limiting behavior of 
the correlation and response functions at large times. These function depend on two different times. When these times 
are large but their separation is small compare to the individual times, the behavior depends only on the separation 
and thus TTI holds as well as the fluctuation dissipation theorem (FDT). On the other hand when the separation of 
times becomes very large TTI and FDT break down although a generalized form of FDT still hold. The breakdown of 
TTI and FDT is one of the main the characteristics of the glassy phase and is referred to in the literature as "aging" . 
Although the dynamical equations are valid for random potential with quite general correlations, in this paper we 
solve the equations for the case of short ranged correlations of the random potential. This case includes the case of 
5-correlated potential i.e. the potential at different points in space are uncorrelated, but we include the case of a short 
correlation length and the case that there is power law decay of the correlation with large enough power. In this case 
the equilibrium solution involved a 1-step RSB as found in Ref. ^2 • We also consider in an Appendix the case of long 
ranged quadratic correlation of the random potential which is exactly soluble without the variational approximation. 
The case of the solution for other long ranged correlated potential will be considered elsewhere. For the statics this 
case involves continuous RSB |6(, and thus the dynamical ansatz should be different. 



The model is defined as follows. The Langevin dynamics is assumed to be governed by the Hamiltonian if[x], 



II. THE MODEL 




(1) 
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where x(s, t) is a d-dimensional vector representing the position of chain bead s at time t. Beads are numbered 
continuously from s = to s = L. rj{s, t) is Gaussian noise: 



(r ] {s,t)r 1 {s',t')) T = 28{s-s')8{t-t')T 



(2) 



due to contact with a heat bath at temperature T. This dynamics is the simplest dynamics for a polymer chain, 
referred to in the literature as the Rouse model 22] . The Hamiltonian H[x] = Hq[x]+ H ran d[x\ contains a deterministic 
part Hq [x] and a random part H ran d [x] . The Hq [x] is defined as 



H„\x] = I 



ds[M(dx(s, t)/ds) 2 + fj,x(s, tf 



(3) 



where M — d T/b 2 K , d is the number of spatial dimensions and bj< is the Kuhn bond length of the polymer. This 
representation is the simplest representation of a polymer as a gaussian chain in the continuum approximation. The 
parameter /i plays the role of a finite volume since the polymer is confined by the harmonic potential to a finite 
region of space. Thus fi — > is the large volume limit and | In /x| oc In V for a volume V. |12| The random part H ran d 
describes the interaction between each bead and the external random potential: 



- f 1 

JO 



dsV(x(s,t)). 



V(x) is a short-range potential, and for simplicity we take it to have a Gaussian form, 

{x-x'f^ 



(V(x)V(x')) 



(dna 2 ) d / 2 



exp 



da 2 



(4) 



(5) 



d is the dimensionality of the system, and a parameterizes the range of the potential. In particular, for a ~ > 0, 
(y {x)V (x 1 )) — > gd(x — x'), and we recover the potential used in More generally we can take 



(V{x)V(x')) = -dJ 

for some function J(z). For the case represented by Eq.JSJ), 

J(z) 



l\2 



d{dira 2 ) d /' 



(x ~ x') 



■ cxp(— z/a 2 ). 



(6) 



(7) 



This model admits a stationary solution characterized by a Gibbs distribution. The equilibrium partition function 
for this solution is given by 



Dxe 



~2T Jo ds[M(dx(s)/ds) 2 + ux(s) 2 + V(x(s)] 



(8) 



III. MAPPING TO THE FIELD THEORY 

Here we will follow closely the notation of Ref. [2C|. Using the standard Martin-Siggia-Rose formalism j^], the 
dynamical average of any observable can be calculated as 



(0[x,x]) T = I DxDxD^D^O(x,x)e' 
with the following dynamical action: 



-S[x,$4,(] 



S[x,x,£,£,] 



dtds 



,d,, dH\x] 
-Tx{s, tf + x{s, t) ( — x(s, t) + a _,_ l ^ 



dx(s, t) t 

- J dtdstts, t)lt(s, t) + J dtdsds'Us, t) dx **^ t) W, t) 



(9) 



(10) 



i, £, £ are auxiliary fields which appear in the formalism. To make for a more compact notation we introduce the 
superficld <£>: 



$(s, h, 0i, 6i) = x(s, h) + e(s, ii)0i + 6^(s, ti) + heix{s, t t ), 



(11) 
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where 8 and 8 are Grassmann variables (anti-commuting c-numbers) . For X, X' e {8,8,8' ,8'}, {X, X'} = and 
/ dXX — 1, the rest of the integrals being zero. In the following, for practical reasons, the more compact notation 
$(s, 1) = $(s, t±, 9i, 6i) will be used. Also, the integral symbol J d8\d8\dl\ will be denoted by J dl. 
In supersymmetric (SUSY) notation Eqs. J3J and (|10|l translate into (|12[1 and 

(0[*]) r - /" D$0[$]e- s[ * ] , (12) 
=5 [$]+5 rand [$], (13) 

where 

So[$] = (1/2) /" dsdlds'd2$(s,l)ii: l s 2 s '$(s'2), (14) 
5 ra „ d [$] = y dldaV(${8,l)), (15) 

and 

i^f ' = S 12 Ss S 'K! , K{=T [n/T - (d/ds) 2 ] - , (16) 
081081 o8\dt\ ot\ 

As noticed by De Dominicis |3^| the expression in Ea. i|12|) is already normalized, so the average over the quenched 
random interactions V can be done directly on i|12|) : 

((A[$])t)v = J D<i>A[^]e- (So[ ' s ' ]+Sl[ ' s ' ] \ (18) 

where exp(— Si[$]) = (exp(— S ran d[&]))v- The average over V can be done easily, leading to 

*[*] = \ Jdsd S 'dld2J ^(^l)-^,2)) a ) . (19) 

The dynamical action Sav = So+Si closely resembles the effective Hamiltonian obtained in the static replica approach 
of refs. This rather general similarity between replica and SUSY treatments has been discussed in ref. |29| . 

Instead of summation over replica indices in we have J dl J d2. 

IV. VARIATIONAL APPROXIMATION 

Since the model cannot be solved exactly, we proceed by using a variational approximation, first introduced by 
Feynman. We Assume that fields <I> are approximately described by a quadratic action 

S var = - y dldsd2ds'$(s, l)G(s, 1; s', 2)- 1 $(s', 2). (20) 

This approach has been widely used in statics. Here we apply it to a dynamic calculation. The goal is to calculate 
Fdyn , given formally by 

e -F dyn =e -({S AV -S var )) vare - iva% (21) 

where 



e 



-J W = J me -s va .r t { )var = e f w J £,$(.) e -S w . (22) 



In usual statics, for problems without disorder, the variational approach is related to a maximum principle. The 
equivalent of Eq. lf?T|) leads to the inequality 

e -F > e -((S AV -S var )) vare -f var _ (23) 
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In the present dynamical problem, as well as in the static problem with replicas, such a maximum principle is not 
known, and the variational free-energy cannot be claimed to be an upper bound to the true one. Despite that, the 
variational approach has been argued to give exact results for directed manifolds in the limit of infinite embedding 
dimensions [34. l35|. giving some justification for its use even at finite dimensions. For real polymers we obtain a 
meaningful solution for 1 < d < 4 and we cannot use the large d limit directly. 

The dynamical variational free-energy Fd yn — ((Sav — S var ))var + fvar is given by 

TP, — i p( 2 ) i p( 3 ) (OA) 

Fdyn — r dyn T £ dyn T r dyn> 

Fdyn = ~/ dsdlds'd2K{iGli (25) 
*fi = -^rlnG (26) 
F^Ud^/jf^l^^)) . (27) 



with 



(3) 

We proceed to calculate the last term F\ ' . Using the identity 



< J(($ - $>'?/d)) var = J d d yJ{y 2 /d) J -0^ exp(-zp • y) (exp( ip • (* - <&')}_ , (28) 



it is easy to verify that 



(exp( ip • ($ - *')) wr = exp(-lp 2 i?-'), (29) 



where 

B{i = G(s, 1; 8 , 1) + G(s', 2; s', 2) - 2G(s, 1; s', 2). (30) 

Defining 

J(a) = J d d y J(y 2 /d) J expHp • y) exp (~sfj 

= T^L°°dx* d/2 - 1 e- x J (31) 
we observe, by substituting Eq. lt2*9"jl in Eq. (|2*gjl that 

(j(($-$') 2 /rf)) Mr = i(<)- (32) 

Thus 

= I / dsdld^J^l'). (33) 



For the case that J is given by Eq.JTJ) we find 



da 



2 \ -d/2 



d(27r) d / 2 V 2 

and recall that a — > for a <5-correlated potential. If J(a) is of the form 

2(1-7) 

for large a, i.e. it involves power law correlations of the disorder, then 



J( a ) = 7^ 7, (35) 



J(a) = (36) 

for large a, where [34[ 

7 = 7 if 7 < 1 + d/2 (37) 

7 = 1 + d/2 if 7 > 1 + d/2 (38) 

g (x g. (39) 
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V. EQUATIONS OF MOTION IN SUPERSYMMETRIC NOTATION 



Given the Fd yn , one can derive the equations of motion from the stationarity condition 
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5G% 



The most complicated term is ^ , F^ n . From it is 



<5g;;' d v r 



- I d3d4dudvJ ' [B%%] {S SS >S US S 13 S 23 + S SS >S VS S 14: S 2 4 - 5 US S VS ,S 13 S 24 - S US >S VS S 14 S 23 ). 



Eq.lgl} simplifies to 



6 -F (3) - d 

=' r dim ~ U 



§G ss> dyr, 



5 SS ,5 12 / dUuJ'iBlD-J'iBH ) 



The variations of F$ n and Fr' n are trivial. Using Eq. JUJJ and (EH leads to 



*3 ' -(Gf 2 ')- 1 +2 



5 a8 /5i 2 / dMuJ'{Bll)-J'{B{l 



0, 



which can be written as 



K{G\l = S 12 5 ss , + 2 j d3duJ'(Bl1)(G^ - G\l ). 
Due to translational invariance in the variable s, G depends only on the difference s — s' . Thus 



GSS /~1S — S 

12 — u l2 



ft is useful to define following Fourier transforms 



Gi 2 = 



— iks /~ik 



Gk 
12 



Since < s < L the corresponding wave numbers k take the values k — (2tt / L)n where n — 0,±1,±2, • 
following it will become necessary to separate the k — component from k ^ 0. Also 



G~ k 
19 



dse iks G% 



Then Eq.|@2} translates into 



fj, + Tk z - D\ 



(2) 



G k = 6, 



where 



12 =o 12 ^2j d3duj'(B^)(e lku G k 32 -G k 12 ), 
^ 3 = 7E^n+^33-2e- lfc ' M Gj 



T 13) 



(40) 

(41) 
(42) 

(43) 
(44) 

(45) 

(46) 
• •. In the 

(47) 

(48) 
(49) 



VI. DISENTANGLING SUPERSYMMETRY 

G\ 2 encodes 16 correlation functions, out of which only two, correlation and response function, are independent 
and nonzero: 



((x(s, tx)x{a\ t 2 )))/d = C(s, h; s', t 2 ) = ~ £ e tt <—'>C fc (ti, t 2 ) 

k 

{{x(s, taxis', t 2 )))/d = R(s, h;s', t 2 ) = jY, e^ s - s '^R k {h,t 2 )- 



(50) 
(51) 
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Also, by adding an external field term to the original Hamiltonian H[x] — > H[x] + J dsdtx(s, t)h(s,t) one gets 

«x(Mi)z(s',i2)» = 6h{ s>,t 2 ) {{xis > h))) ' (52) 
i.e. R(s, t\\ s', t 2 ) describes the response to an infinitesimal field applied at time t 2 and bead s' . Thus, G k 2 reduces to 

G k 2 = C k {t u t 2 ) + (0i - 6 2 ) [OiRkfah) - e 2 Rk(ti,t 2 )} . (53) 

It follows that 

B\ 2 = B u {hM) -\Y, e ~ tk ' U ^ - ^) {°i R k'(h,h) - 6 2 R u ,{t u t 2 )] , (54) 

k 

with 

B u (t 1 ,t 2 ) = (((x(u,t 1 )-x(0,t 2 )) 2 ))/d= [Ck(ti,t 1 ) + C k (t 2 ,t 2 )-2e-* ku C k (t u t 2 )] . (55) 

k 

After disentangling the equations of motion in SUSY notation (see Eq. I48fl by using (|53I55|) gives 

\pL + Mk? +d/dt]C k {t,t') = 2TR k (t',t)+2 f dt a f duJ'[B u (t,t 3 )]R k (t',t 3 )e tku + 

Jo Jo 

+4 f dt 3 f duJ" [B u (t, t 3 )] R u (t, t 3 ) [C k (t, t') - e* ku C k (t 3 , t')} , (56) 



jo 



t /-L 

k'b 



[,j + Mk 2 + 8/dt]R k (t, f ) = 5(t -t')+A dt 3 / dud" [B u {t, t s )] R u {t, t 3 ) [R k (t, t') - e lku R k (t 3 , t')] , (57) 

Jo Jo 



where we defined 

R u (t,t') = jY, e ~ lkUR ^ t > t 'y ( 58 ) 

k 

VII. ANSATZ FOR THE CORRELATION AND RESPONSE FUNCTIONS 

These equations of motion are coupled integro-differential equations which in principle can be solved; the initial 
conditions are given by C' k (0, 0) and we use Ito's convention R(t + e, t) — > 1 as e — > from above. It is well known 
that asymptotic solutions of such equations can be characterized by few parameters and it is possible to solve those 
equations analytically. El El E III III El 

For t, t' -> 00, r/t' < 1 and r = t - t', TTI holds 

Urn C k {t,t) = q k , (59) 

t— »oo 

lim C k (t + r,t) =C k {r), (60) 

t — >oo 

lim Cfc(r) = q k , (61) 

T — >00 

and 

lim R k (t + T,t) = R k (r). (62) 

t — >oo 

We will refer to this regime as the stationary or TTI regime. In addition to the TTI regime, there is another long 
time non trivial regime, characterized by t, t' — > 00, fixing a = h(t')/h(t) and < a < 1, where the function h(t) is an 
increasing function of t which the asymptotic analysis performed here is not able to determine. In this aging regime 
one has 

lim C k {t,hr\ah{t))) = q k C k (a), (63) 

t — >-oo 

lim q k C k (a) = g ,fc, (64) 
lim 4(a) = 1, (65) 
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and 



lim R k {t,od) = -Rk(a). 

t — >oc t 

Also, for future convenience, it is useful to introduce the following order parameters: 



(66) 



~b(u) 
b(u) 
b (u) 



Ed 



— iku\ ~ 



k — e qa.k) 



(67) 
(68) 
(69) 



Also 



iT(a) 



^fe- e - lte C fc (r)], 
^fe- e - lfc > fe C fe (a)], 



-iku 



R k (a). 



(70) 
(71) 
(72) 



VIII. EQUATIONS RELATING ASYMPTOTIC VALUES OF CORRELATION AND RESPONSE 

FUNCTIONS 

Using the ansatz discussed in section IVIII one can derive the following equations for C k {t,t') in the TTI regime: 
(fi + Tk 2 + d/dT)C k (r) = 

2TR k (-T) + | f duJ'[b(u)\ [C k (r) - e lku q k ] 
1 Jo 

-1 [ duJ'[b(u)](l - e^)C k {r) - 1 [ dr' [ due^J'[B^r - r')}^0± 



o ± Jo Jo 

+2 f dp f L due tku J'[B u (p)}R k (p)+4: f dp [ ' duJ"[B u (p)]R u (p) \c k (r) - e lku q k Ck{p) 
Jo Jo Jo Jo L 



(73) 



It is also possible to derive similar equations for Rk(r) which, due to the Fluctuation Dissipation Theorem (FDT) 

1 dC k (r) 



Rk(r) 



T dr 



(74) 



are completely equivalent to Eq. J73J). 

In the aging regime one gets the following equation for qkC{a): 



[p + Mk 2 -^{ dpi due zku J"(B u (p))R u (p) 
Jo Jo 



■1 r-L 

q k C k {a) = 2 1 dp due lku J'{B u (ap))R k {p) 
'o Jo 



+ / due iku J'(B u (a))(q k -q k )-A dp due* ku J"(B u (p))R u (p)q k Ck(p/a) 
1 Jo Jo Jo 

pi pL poo pL 

-4/ dp du^ ku j"{B u {p))R u {p)q k Ck{a/p) + A dr' du J"[B u (t')]R u (t')(1 - e lku )q k C k {a). (75) 

Ja Jo Jo Jo 
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For R k (a) we obtain, 



p + Mk^-i dp duJ"{B u {p))R u (p) 



1 dp 



Rk(a) = - 



T 



due lku J"(B(a))R u (a)(q k -q k ) 



-4 I — I due lku J"(B u {p))R u {p)R k {a/p)+4 dr'\ du{l - e lku )J / '[B u (T , )}R u (T , )R k (a). (76) 

Ja P Jo JO JO 

Again, one can see that both Eq. i(75() and Eq. 176(1 can be solved by the ansatz 

r> f \ x c dC k (a) 

Rk(a) = Tflk—, ■ (77 

1 da 

Eq. 177(1 is commonly referred to as a generalized FDT (GFDT). The parameter x c corresponds to the corresponding 
parameter in the static replica solution. In the context of replicas it was first introduced by Parisi and should not be 
confused with a spatial coordinate. This parameter must satisfy the inequality x c < 1. In principle, Eq. ((77|) could 
have been written as 



a , v x ck (q k C k (a)) dC k (a) 
Rk{a) = q k - 



T 



do 



(78) 



which could be applied to a many-step RSB scheme. However, as for the case of directed polymer with short ranged 
correlated random potential we found jl2j that a solution with one step RSB is appropriate, and it is sufficient to use 
the simpler ansatz given in Eq. ((77(1 . 
For t = t' and t — ► oo Eq.JSSJ) gives 



T 



(fi + Mk 2 )q k = T + — duJ'[b(u)](q k - e iku q k ) - - / du(l - e lku )J'[b{u)]q k 



T 



+2 dp due iku J'[B u (p)]R k (p)+A \ dp j duJ"[B u (p)}R u (p) [q k - e iku q k C k (p)] 
Jo Jo 

Eq. 1(73(1 for t — > oo and then r — > oo results in 



(79) 



(p + Mk 2 )q k = 



T 



du{l-e lku ){J'[b{u)]-J'[b{u)\}qk 



+ - j due iku J'[b(u)](q k -q k ) + 2 j dp j due J 1 (B u {p))R k (p) 



dp / duJ"(B u (p))R u (p)q k [1 - e lku C k {p)} 



Also, Ea. (|75|l for a — > gives 



(80) 



(p + Mk')q . k = 2 / due tku J'[b (u)] / d P R k (p) + 



TJo 



due lku J'[b (u)](q k -q k ) 



2 
T 



du(l - e lku ){J'[b(u)} - J'[b(u)}}q , k + 



T 



du(l-e lku )J"[B u (p)]R u (p)q „ 



(81) 



Eqs. JZHl, and JSU and their origin Eqs. JHJ and contain both TTI and aging parts. Thus, in principle, 
there are two anstze for solving them, leading to two physical behaviors an ergodic phase (this name was coined by 
KHF [2(| in the sense of non-glassy) characterized by TTI and FDT where the aging behavior is completely missing, 
and a glassy phase containing both stationary and aging behaviors. 



IX. ERGODIC PHASE 



By An ergodic phase we mean that the external parameters are such that only the stationary solution exists and 
not the aging solution. This happens when the only solution is with q k — qo tk . Technically, assuming that aging is 
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absent amounts to setting R k (a) — and O k (a) = 1 in (|79ll , 1801) and i|81|l . (Equivalently, one could start from 15611 
and l|57l) and exclude the aging part from the beginning, leading to the same equations.) Thus, in the ergodic phase, 
equations {75). (|%0f and (jHU reduce to 

(^ + Mk 2 )q k =T+^ [ duJ'lb(u)}(q k ~e lku q k )-?- [ du{l - e lku ) J'[b(u)]q k , (82) 
1 Jo 1 Jo 

(fl + Mk 2 )q k = %\ L due lku J'[b{u)]{q k - q k ) + | C du(l - e lku ){J'[b{u)] - J'[b(u)]}q k , (83) 
1 Jo 1 Jo 

(fx + Mk 2 )q Q . k = | ^ due tku J'[b (u)](q k -q k ) + ^ £ du(l - e tku ){J'[b(u)} - J'[b(u)}}q , k . (84) 

Note that l|83|) and l|84|l enforce q k = qo yk which is just equivalent to C k (a) = 1, so one gets only two equations. 
In order to solve these equation we define the following constants, which are themselves functions of q k and q k : 

E k = %\ due lku J'[b(u)}, (85) 
1 Jo 

Fk = ^l du{l-e^)J'[b(u)]. (86) 

In terms of these constants the equations for q k and q k become: 

(M + Mk 2 )q k = T + E a q k - E k q k - F k q k , (87) 
(M + Mk 2 )q k = E k (q k - q k ) + {E - E k - F k )q k . (88) 

The solution of these equations is 

- _ Tjji + Mk 2 -E + 2E k + F k ) 

qk (n + Mk 2 -E + E k + F k ) 2 ' l8yj 

Qk = {n + MV-Eo + Ek + Ftf- (90) 

We now show that the ansatz E k = for k ^ solves the equations. Since F = we find 

T TE , . 

qo = ~ + (91) 



T 

fi + Mk~E + F k 
TE 



- - , F _ , E ( 92 ) 



qo = — —, Qk^o = 0. (93) 
A* 



Using these solutions we see that 



is independent of u. Thus 



IT IT 1 



2 r 2 r 

£fc = J^(6) / ^e lte = — J'(6)4,o- (95) 



which validates our ansatz. Thus we can write 



2L f/ / 2T 2T v-^ 1 | 



Also 

9 9T 1 — p—iku 

£ ^ L ^ fi + Mk 1 - E + F k 
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and hence 

* -| [ *(. - e-)i' E „ + ^" +JV ) <<*> 

This equation together with Ea. (|96|l gives a complete set of equations in the ergodic case. They are identical to 
Eq.(4.15) and (4.16) in Ref. |6| derived by the replica method for a quantum particle in a random potential(the 
notation there is slightly different but it is easy to identify the corresponding variables). Thus in the ergodic case 
there is a complete agreement between the dynamical calculation and the replica calculation. 
From the results above we obtain 

q = lim ((x(s,t)x(s,t)))/d = ~ V g fe = — + —j + — V — • ^ ~ 

t^oo L k L f 1 L r L J fi + Mk z - E a + F k 

1 T E 

q = lim \im((x(s,t + T)x(s,t)))/d= - V % = — -2. (99) 



X. SPIN GLASS PHASE 

In this phase go ^ go,fc and there is a time regime where the aging behavior takes place. This phase corresponds to 
the RSB phase of the statics. We introduce the functions 

E ,k = y [ due iku J'[b (u)], (100) 
1 Jo 

Keeping the aging parts and using the GFDT, Eqs. I|79(l . I|80JI and (|81|l can be transformed into 

(p + Mk 2 )q k = T + E (l - x c )q k - E k (l - x c )q k - F k q k + E ,oX c q k ~ E , k x c q 0>k , (101) 
{H + Mk 2 )q k = E k (q k - q k ) + [25b (1 - x c ) - E k (l - x c ) + E Q:0 x c - F k ]q k - E 0tk x c q 0tk , (102) 
(/i + Mk 2 )q 0tk = E 0ik q k - E 0)k (l - x c )q k + [25 (1 - x c ) - E k (l - x c ) + E , x c - 2E 0tk x c - F k ]q , k . (103) 

In this case, similar to the ergodic case we will use the ansatz E k = and Eg }k = for k ^ 0. We will see that this 
provides again a consistent solution. Using this ansatz the solution of Eqs. I|101|) . I|102|l and [jl03|l for q k , q k and go,fe 
becomes 



T 

" x /i + MAr - 25 + S + F fc 

T(^25 + j5 , S) 

?o = 2 ^ T^i ' = 0, (105) 



where 



Using this solution we see that 



Qo.o = TE y'° qo,k^o = 0. (106) 



£ = a; c (25o - 25 ,o)- (107) 



2T 2T 1 

= Zo7T^) + - g ^W^TETn 55 6 ' (108) 

+ S) L £f o /j, + Mk 2 - E Q + Z + F k 



both independent of u. Thus 



2 f IT 

E k = ^J'(b) J due lku = — J'(b)S k<0 , (110) 

2 ilfi. \ f L J..„iku fl 



E ,k = ^J'(fco) / due tKU = —J'(b )d kfi , (111) 
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consistent with our ansatz. If we denote 



-ik u 



a{u) = T £ M + A/fc2-^ + E + F fc ' 



we can write 



E = p' 



o(0) 



2T 



o(0) 



2T 



1 



and thus 



E = 



*2Lx c 
T 



a(0) 



2T 



J' 



i(0) ■ 



E 
2T 



E 

flX c 



For Ffc we get instead of Eq. (JSSJ 



F k = - I du(l-e tku )J'[a(0)-a(u)} 



Let us dehne the functions 

[ L due lku J"[b(u)]R u (l), 
1 Jo 

from Ea. (|77|l and the fact that qk^o = it follows that R k ^o{a) = 0, and hence also Dk=£o — 0. Thus 

Dk = |j"[b]i?o(l)4,o. 
Furthermore, Eq. (|76|l with a = 1 and fc = gives 

+ S)-Ro(l) = -2D {q - qo) 

which can be written as 



/i + E 



J"(b)Ro(l), 



= R Q (1) 



1 



Eq. i|120|) with J2q(1) + 1 implies the marginal stability condition 

2T 



J"(6) 



l(a» + S) 



J" 



o(0) 



2T 



L(/i + £) 



(112) 

(113) 
(114) 

(115) 

(116) 
(117) 

(118) 
(119) 
(120) 

(121) 



Eqs. I|115|) . (|144l) . (|145|l and l|121|) fully solve the model. These equations, with the exception of Eq. ll21|) . are the 

same as the equations for the statics that we obtained using the replica method in Ref. 6] and Ref. |12(. Equation 
(|121J) though, is of pure dynamic origin, and differs from the corresponding equation obtained in the replica method 
by extremizing the variational free energy with respect to the RSB breakpoint x c . We give here for comparison the 
corresponding equation obtained in the statics replica calculation |12| : 



T 2 2 
T 2 



J 



O(0) 



2T 



2La 



T ^(/t + £ 



•J' 



a(0) 



L(fi + E 
2T 



1 



- J 

E 



o(0) 



2T 



£(A» + E) 
E 



/i + E 



In I 1 + - ) = 



(122) 



Equation l|121|) can be obtained in the replica calculation by requiring marginal stability i.e. the condition of a 
vanishing replicon mass, but it does not correspond to the optimized variational solution. 
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XI. SOLVING THE EQUATIONS 



In this section we discuss the solution to the equations derived in the last section. Although a full solution can be 
found numerically, our goal here is to go as far as one can with analytic methods because it gives a better understanding 
of the nature of the solution. An analytical solution becomes possible for a system of large volume (small value of li) 
and whwn the polymer is very long (large L). 

Before we consider the solution of the equations derived above for a polymer, let us first discuss the special limit 
where the problem reduces to a classical particle in a random potential. Looking back at the original hamiltonian one 
realizes that in the limit M — > oo and T — > TL (or alternatively L = 1) the problem should become the same as the 
particle problem discussed in Refs. 0, |3|| |3tJ ■ Looking back at our equations we see that in the limit M — > oo only 
the k — component of qk survive. 

In the ergodic case Eq is given by 



E = -J'[ — ), (123) 



and the expressions l|99|l for q and q become 



T V [i ) ' 



9=- + ^, (124) 

a = (125) 



These are exactly the same as equations (57)- (59) derived by Engel .36} using the replica method and correspond to 
his replica symmetric solution. 

In the spin glass case we again see that as M — * oo we have a(u) = 0. Thus Eas. (|113|) . H114fl and l|115l) with a(0) = 
and L = 1 agree exactly with equations (74)-(76) of Engel [3|J using replicas with 1-step RSB . However it is Ea. (|122|) 
with a(0) = and L = 1 that agrees with Engel's Eq.(78) and not the dynamical equation (|121H with a(0) = that 
agrees with the dynamical equation derived in Refs. |18l l37l |. We also find that 

T f l+ E 1 + E^\ (126) 



agree with Engel's Eq.(79). 

Let us review Engel's solution to the variational equations for the short range, random potential case in the limit 
of small fx (but still fi ^ 0). Ea. (|122|l with a(0) = and L = 1 gives 

xl t {2T 

Thus 



r2 -Jl — ) -ln( M )=0, (127) 



, c = Tfi^V /2 , (128) 

where we used the fact that 2T/S — > for [i — > as will become clear shortly. Notice that in this case of a particle in a 
random potential we must consider a random potential that is regularized at small distances (not a strict S- function), 
so J(0) is properly defined. Equation 1115|l then gives 



2x c Pr^ n t//^ ( I ln (M) I \ 



1/2 



Also E'o o vanishes like 



(\ (2+d)/4 
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and E = 2J'(0)/T. Thus 

. TE |J(0)|V 2 



(131) 



The importance of Ea. l|131[l is that it is consistent with the following Imry-Ma type argument: Since the random 
potential has a gaussian distribution of variance g, in a volume of radius r g the minimum of the potential is given by 
V m ~ — yj g \n(r g ). This result follows from the fact that when we pick r g numbers from a gaussian distribution with 
zero mean and unit variance, the lowest expected value of the potential is given by the solution of the equation 



/' 

J — ( 



1 

exp(-V 2 /2g)~-. (132) 



Thus the Hamiltonian is H(x g ) ~ (l/2)/ir 2 — yj g In r g and minimization gives 



VM|ln/i| 1/4 ' 



(133) 



yielding 



g=((^ 2 ))- 9 '[ 2 i ( 134 ) 

This is a nonanalytic expression that cannot be obtained from perturbation theory. The critical temperature below 
which the RSB solution is the stable one is given by the condition x c = 1 or T c = (J(0)/| ln/iQ 1 / 2 . 

In the dynamics approach though, Eq. (|121f) replaces Eq. (|122f> . For small since E 3> /i Eq. (|121|l gives 

- 1 = ^"(0), (135) 

or 

E = 2|J"(0)| 1/2 . (136) 
again Eq = 2J'(0)/T and _E .o vanishes like ^( 2+d )/ 2 and thus from Eq. (|115fl we get 

. JJ"(fl)\ 1/2 



J'(0) 

We also find that 

- TE n J'(0) 



m I J" (0)|V2 



(137) 



(138) 



We observe that in the dynamical formulation ln/x is replaced by the constant J"(0)| J(0)|/(J'(0)) 2 which is equal to 
(d+2)/d for the short ranged correlated potential. The Imry-Ma result is not satisfied, even though the distance of the 
particle from the origin still diverges for fi — ► 0. This is probably due to the fact that in the dynamical formulation 
the particle is impeded by large barriers to search the entire volume as effectively for the lowest minimum of the 
random potential as is obtained in the statics, even in the limit of large times. The fact that the dynamical solution 
differs from the statics in this case was observed in Ref. but they did not compare the solutions in detail. They 
observe that in the dynamical solution the free energy is higher than the replica free energy obtained in the statics. 
This discrepancy may also be due to the fact that the initial conditions of the dynamics are completely random and 
are not weighted with an appropriate Boltzmann factor . 

We now return to the polymer problem and discuss the solution to the previously derived equations for the case 
of a random potential with short range correlations characterized by7 = l + <i/2. The ergodic phase for large times 
in principle should correspond to the replica symmetric solution of the statics. In the equilibrium solution such a 
solution arises in two situations: First, in the infinite volume limit, when /i = 0. In this case the polymer totally 
collapses to a size of a sing le monomer, since deep minima of unbounded bottom exist in this limit. This case was 
discussed before in Ref. |12J and it was shown that the solution corresponds to the case of an annealed disorder. Since 
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the solution was discussed before we will not repeat it here. In this case the equations for the dynamics are the same 
as obtained in the replica solution. Second, for the finite volume case the ergodic solution is applicable when L < L c 
where L c is the value for which the spin glass solution ceases to exist, see below. It turns out that L c is a number of 
order unity, thus except for very short chains the non-ergodic solution is the correct one. It is a subtle question if the 
domain of validity of the ergodic phase in the dynamics and in the statics coincide. For the case of infinite volume 
one has to take the limit /i — > before the limit of large time, we will not discuss this question further here. 

In the Appendix we give the solution of the ergodic equations for a special case of a long range random potential with 
quadratic correlations and show that it gives the same solution as previously obtained by Shiferaw and Goldschmidt 
|13| using a different method. Of course for a quadratically correlated potential the variational approximation is exact. 

We proceed to discuss the solution to the equations in the spin-glass case when the random potential has short 
range correlations (7 = 1 + d/2). As mentioned above, when the volume is large but finite the spin-glass solution is 
the appropriate solution when L is large enough [l2]. For the spin glass case we have seen that one of the equations 
of the dynamics differs from one of the equations derived in the statics using the replica method. We now re-derive 
the solution of the equations of the statics and then show how the solution of the equations of the dynamics differs 
from it. The reason that we reconsider the solution of the equations of the statics is first because the variational 
scheme that emerged above is more general than the scheme employed in Ref. |l2j since it involves more variational 
parameters like the scheme used in Ref. [(J, and we want to show that the end results are the same. Second, the steps 
of the solution will facilitate the solutions of the dynamical equations. 

In the limit of small fi (but still fj, ^ 0). Ea. (|122[) gives 



where we used the fact that 2T/(L£) <C a(0) for large L ( and even for finite L as n — > as can be checked a posteriori 
with full solution). Equation 1|115|) then gives 




(139) 



Thus 




(140) 




(141) 



and from Ea. (|114|l one obtains 




(142) 



It is convenient to define the variables 



A/, 



E +'E + f j,, 



(143) 



Using Ea. l|113|) . A& satisfy the equation 




(144) 



and in terms of A& the constant a(u) is given by 




(145) 



where the last expression is valid for large L. The integral on the right hand side of Eq. I|144|) converges as L becomes 
large and thus to leading order the parameters Xk are 0(1) with respect to L. To simplify the integral further notice 
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first that the integrand is invariant under the transformation u — > L — u thus it symmetric about u — L/2. The total 
integral is twice its value up to u — L/2 and A& can consistently be taken as real. For large L we obtain 



A fc = p + S - Jj"[a(Q)] + - / du (1 - cos fcu) { J' [a(0) - o(«)] - J' [a(Q)]J + 0(l/i), 



where £ is given by Ea. ll41|l and 



a(u) = 2T 



dk cos ku 



, 2ir Mk 2 + X k 

Note that we have added to the integral in Ea. (|146(l a term 

4 9n 



T 



J>(0)] / dwcos(fcu)=0, k ^ 0. 



(146) 



(147) 



(148) 



Our goal is to characterize the behavior of the end-to-end distance of the polymer that can be extracted from the 
correlation function 



b(L) = lim (([x(L, t) - x(0, t)} 2 )) = 7 Y> - e- 4fci )% = o(0) - a(L) 



(149) 



fe^O 



We are now going to argue that as /i — > (but still finite), Afc satisfies, to leading order in | ln^t| the scaling form 



X k = (, 9 | lnp\) S X + y^f-f (k(9\ ln^|)- 5/2 ) , k > 



(150) 



where 5 = 4/(4 — d) and X- k = Xk- For small k the function / satisfies f(x) ~ a; 2 and is regular around 0. For large 
x, it can be shown that f(x) ~ x d / 2 . Substituting Ea. (|150(l into Ea. (|147(l and changing the integration variable 



we find 



a(u) 



2T 



k^k(g\lnp\) 5 / 2 , 

dk cos(k(g\lnp\) 5 / 2 u) 



(g\ \np\y/ 2 2ir Mk 2 + A + jj^fik) 
From which it follows that to leading order in | hxfi\ 



a(u) 



T 



AM 



= ( 5 |ln M |)- 5/2 e-VA/M( s |in M |)^^ 



(151) 



(152) 



(153) 



and in particular a(0) = T/y MX (g\ ln/i|) & l 2 . Substituting on the rhs of Ea. ll46|l and changing the integration 
variable 



u^u(g\lnp\)- S ^(M/X)^ 2 , 



(154) 



we find for small fi 



X k = 



(2n) d / i 



T 



-(d+4)/4 



MX. 



(^"" 1W4 /(*(^,D^ 



In /'I 



where 



2M 



T \ 



-d/2 



(2tt)<V 2 T 2 \y/Mx) 



du 



{l-cos[,(M/A)^]}[-_i^ 



(155) 



(156) 
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we see that for consistency 5 must satisfy S = 4/(4 — d). We also see that the parameter A is given by 

4/(4-d) / , ,-\ (4+rf)/(4-rf) 



(27r) d / 2 ; V T V 



(157) 



Using Ea. l|153|) we have 

b(L) = -^=(pjln/x|)- 5 / 2 fi-e-VWfel 1 ^!)^ 2 ^ . (1 58) 
V AM v 7 

Thus the end-to-end distance is given by 

i? 2 ~-J=(.9|hiMir 5/2 , (159) 
V AM 

from which it follows that 

R~ ( 5 |ln//|)- 1/(4 - d) - ( 5 lnV)- 1/(4 - d) . (160) 

Here V is the total volume, and this result has the correct q dependence and the subtle InV dependence as was 
originally argued by Cates and Ball [TlJ and derived in Ref. [13 using the replica method. 

The important observation is that when k increases from to \J\/M(g \ In /i)^ 2 , Xk changes only by a factor 
1 + 0(1/| In fi\) and thus the results of the single A parameter used in the variational scheme of Ref. [12J remain intact. 
The effective mass of the low lying non-zero modes is approximately given by \fk~k ~ 

The above scaling arguments are valid provided a(0) > da 2 /2 (see Ea. (|34|) which is what we are going to assume. 
Of course for a delta function correlated random potential a — ► and this condition always holds. For this condition 
to apply in the case of other short range correlated random potentials, g the variance of the disorder, may not be too 
large, such that the size of the polymer is not smaller than the correlation length of the disorder. 

The parameters S and x c are given to leading order in L and | In fj,\ by 

S = rf 4 /( 4 - d )(2^)- 2d /( 4 - d )(M/r 2 )( 4+d )/( 4 - d )(. g |ln A1 |) 4 /( 4 -' 1 ), (161) 



jd-2 



-l/(4-d) 



(162) 



We see that x c < 1 for large enough L, and for d > 2 for fixed L and small enough /1. The value of L corresponding 
to x c = 1 is denoted by L c . The condition x c < 1 can be written as 



InV < L 



,glnV 



(163) 



which means that the translational entropy is smaller than the binding energy in the typical minimum of the random 
potential, thus implying that the polymer is truly localized. It also follows that x c < 1 is equivalent to the condition 



gL (*-d)/2 

T 2 b d K 



(InV) 



(d-2)/2 



> 1, 



(164) 



up to some unimportant numerical factor. 

We now discuss the solution of the dynamical equations. Eq. I|121|) replaces Eq. (|122|) . For small \i since S ^> /1 
Eq. flUl gives 



1 = —J"[a(0)], 



(165) 



E = 2 



J>(0)] 



1/2 



(166) 
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From Eq. i|115|) we get 

LJ'[a(0)] 

and we see that there is no longer any ln/j dependence as in the replica solution. Again x c < 1 for large enough L. 
Using this equation we see that 

_ 2L ( 2T \ gL ( y 2 \ J"[a(0)]\ \ nm 

The equation for A& again reads 

X k = m + S - -J"[a(0)] + - / du(l-cosfcu){ J' [a(0) - a(u)] - J' [a(0)H + 0(1/L), (169) 

Since in this case there in no ln/j dependence we can repeat the scaling of A& with g as before but not with ln/i. We 
find the g dependence as before, thus 

\ k =g^ A -V\ k , (170) 

where A& are independent of g. We also find that since | In fj\ is absent, the parameter M will be shifted by contributions 
from f(k). We can still argue that Afc are independent of L which is an important result, and that they have the g 
dependence described above. This implies that 

R ~ g- 1 /^) (171) 

independent of L. Similarly we find that S is finite and satisfies E ~ ff 4 /( 4 " d ) and x c - {l/L)g- 2 ^ 4 - d \ This last 
condition implies that the polymer is localized whenever x c < 1 which means approximately gL^ 4 ~ d ^ 2 > 1. 

These results still contain important physics, namely that the size of the polymer is independent of its length and 
also it has the correct dependence on the strength of the disorder, but the more subtle In V dependence resulting from 
a sophisticated Imry-Ma type argument is missing. 

XII. DISCUSSION OF THE DYNAMICAL BEHAVIOR 

Here we make some further comments about the dynamical behavior of the polymer in the spin glass phase. Consider 
the motion of the center of mass of the polymer: 



1 

dL 2 Jo 



B C u{t,t') = {{[x C M(t) - X C M(t')] 2 ))/d = 

L pL 

du I du'(([x(u,t) — x(u, t')] ■ [x(u',t) — x(u',t')])) = 







j [C (t,t) + C (t',t') - 2C (t,t')] . (172) 



Let us denote t' = tyy and r = t — t' . In the TTI regime for large ty/ and r -C tyy, as r increases, Bcm{t) increases 
from to 

2 2T 

b {1] = lim lim B C m{^w + r,t w ) = — (q - q ) = — — — — . (173) 

T— tOOtw— KX) L L[fl + Zjj 

In this regime the FDT also holds. Let us compare this behavior with the behavior of a free chain. For a free chain, 
for large t, t' we have 

IT IT 
B C m(t) = — 1 - e-" T ) ^^oo — , (174) 

and thus when disorder is present S is replacing jj, when yn — > and we observe that 6^1) remains finite as /i —> 0. In 
this regime the polymer becomes trapped inside a local potential minimum and we see that for large r ( but still less 
than t\y) there is no diffusion, even though for very small r we expect diffusive behavior. Thus for large enough r 
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there will be a plateau in theplot of Bcm{t) as a function of r. In this regime the short time estimates of previous 
investigations (see e.g. Ref. |25|)] should be valid. 

However, for large and fixed t\y as r becomes sufficiently large, Bcm will leave the plateau and continue to grow 
above until it reaches the value 

2 IT ( Yj \ 
6<°> = lim Bcm (tw + T,t w ) = -{q - q ,o) = ,, , = U + ■ (175) 

The size of the plateau depends on t w . The larger iy^ the larger the plateau, and the larger value of r required 
for Bcm to increase beyond Thus the polymer does not remain trapped forever but eventually hops to another 
minimum of the potential. Notice also that as fi — > 0, b^ ' — > 2T / (fiLx c ) ~ 2Tg 2 /( 4 ~ d ) / fx. This value is independent of 
L and should represent the typical square of the hopping distance of the polymer among different local minima of the 
potential. The larger the waiting time the deeper the local minimum occupied by the polymer and hence the longer 
it takes it to hop to another minimum. From with this observation we still lack an estimate of the time dependence 
of BcM(t,t') in the vicinity of and b^°\ A more detailed calculation is needed to derive the asymptotic growth 
rate in the different time regimes, and this will be left for future work. 

Another quantity of interest is the single segment (bead or monomer) mean square displacement 

B"=°(i,f) = {([x(s,t) - x( S ,t')] 2 ))/d = 

j J2 [Ck(t, t) + C k (t', t') - 2C k (t, t')] . (176) 

k 

In this case the asymptotic mean square displacement in the TTI regime becomes 

bW=2(q-q)=bW+a(0), (177) 
and the asymptotic value in the aging regime becomes 

6(° e ) = 2(q-q ) = b^+a(Q), (178) 

with M 1 ' and M°) defined above and 

a(0) ~ {TbD^^g- 2 '^. (179) 

The quantity b^J g is dominated by a(0) and thus the size of the mean square displacement of a bead is the same as 
the square of the end-to-end distance. The quantity b^Jg remains essentially the same as b^ for small \i. 

The results of this section i.e. the width of a localized state of the polymer and the average distance squared 
between different localized states is the same as discussed in the interpretation of the 1-step RSB solution in Ref . [ 
section VII, with only the subtle dependence on InV missing. 



XIII. CONCLUSIONS 



In this paper we derived the dynamical equation for a gaussian chain in a short range correlated random potential. 
We used a simple Langevin dynamics and we discovered that there are two possible scenarios at large times: an 
stationary regime where the FDT applies and an aging regime where the FDT breaks down at large time separation. 
In the aging regime FDT can be shown to be replaced by a modified or generalized form commonly referred to as 
GFDT, and involves a parameter x c similar to Parisi's parameter for 1-step RSB. Only 1-step RSB is necessary for 
the case of short range correlations^ In the long range case that was discussed in Ref. Q in the equivalent context of 
a quantum particle in a random potential full RSB applies). 

The stationary regime represents the dynamics of a chain trapped in a local minima of the random potential. 
Eventually after very long time the chain can escape from its pinning and hop to another minimum elsewhere leading 
to history dependence and violation of equilibrium theorems. In the long time limit, for a short ranged correlated 
random potential, the dynamical equations become identical to the equations of the statics as derived from the replica 
method except for one equation that involves x c , which is different from the equation derived in the statics using 
replicas. This is probably due to the fact that starting from random initial conditions the dynamics gets influenced 
by large barriers and does not explore the potential landscape as efficiently as to reproduce the statics even at large 
times. Thus the subtle In V dependence of the statics that emerges from an Imry-Ma type argument even for the case 
of a zero-dimensional object in a random potential is not reproduced by the dynamical equations. Our results are 
based of course on the Gaussian variational approximation but the emerging picture is probably valid. 
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APPENDIX A: RANDOM POTENTIAL WITH QUADRATIC CORRELATIONS 

In this Appendix we discuss the exactly solvable case of a potential with long range quadratic correlations [T^| . 
Since the variational approximation becomes exact for such a potential, our dynamical equations should reproduce 
the solution found using the statics and replica formalism. For this case J(a) = (l/2)ga + const, and so is J(a) 
since 7 = 0. Thus J'(z) — g/2, a constant. (The notation here is slightly different from Ref. |b| where we used 
J(a) = 2cra + const, so g — > 4er.) Only the ergodic case applies in this case since Eq = i?o,o- We find 

Eo = (Al) 



T 

F k = ^(l-S kfi ). (A2) 



We also obtain 



* = f 4 ' + ^W' (A3) 



From these results it follows that 



Qk = ^4,o- (A4) 



9 T 1 9 , T 

^ k 

where the last expression applies for large L, and also 



9 ST^ 9 I A c.\ 



The correlation function b(L) is given by 



9=4- (A6) 



b(L) = -= . (A7) 



For small u this function becomes equal to L/M as in the free case. These results coincide with the results obtained 
in Ref. 0. 

Actually in this case one can write a closed form solution for Eas. (|56|) - (|57[) . The solution is 



C t (i.i') = i±8 t ,n (l - e"" - e-"' + e -<"'+''>) + 

R k (t,t') = 9(t - t') e ^+Mki)(t-t>) (A9) 
Indeed for large times C(t — t') becomes TTI and depends only on the difference t — t', and the FDT holds. 
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